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We consider a mesoscopic four-terminal superconductor/normal metal (S/N) structure in the 
presence of a temperature gradient along the N wire. A thermoemf arises in this system even in 
the absence of the thermoelectric quasiparticle current if the phase difference between the supercon- 
ductors is not zero. We show that the thermoemf is not small in the case of a negligible Josephson 
coupling between two superconductors. It is also shown that the thermoelectric voltage has two 
maxima: one at a low temperature and another at a temperature close to the critical temperature. 
The obtained temperature dependence of the thermoemf describes qualitatively experimental data. 



I. INTRODUCTION 



Transport phenomena in mesoscopic superconductor-normal metal (S-N) structures have attracted a great interest 
in last years (see for example reviews [1,2] and references therein). Due to the proximity effect (PE) such properties 
of the normal metal N as the density-of-states etc are changed, and therefore transport properties of S-N structures 
are changed as well. For example, the resistance of a normal wire a part of which is a segment of a superconducting 
loop is a periodic function of a magnetic flux <!> threading this loop. The period of the oscillations is the magnetic flux 
quantum $o = hc/2e. The amplitude of the resistance oscillations SR m depends on temperature T in a non-monotonic 
way reaching a maximum at T m = ciCl, where c\ is a numerical factor and = / L 2 is the Thouless energy [3,4]. 
A lot of publications are devoted to the study of the electric conductance of S-N mesoscopic structures. Much less 
attention was paid to the study of thermoelectric phenomena in such structures. However the situation has changed 
recently. A number of papers has been published in which the results of both experimental [5,6,8] and theoretical 
[11-15,17] studies are presented. It has been established experimentally that the thermoemf Vth in a S-N structure 
with a superconducting loop also is an oscillating function of the magnetic flux in the loop and the amplitude of 
oscillations depends on temperature in a non-monotonic way [5,6,8]. 

In Refs. [14,15,7] the thermoconductance in S-N mesoscopic structures was calculated and measured. These calcu- 
lations generalized the results for the thermoconductivity in superconductors in the intermediate state obtained by 
Andreev a long time ago [16]. The thermoelectric effects in superconductors caused by the thermoelectric component 
jth = i?VT in the quasiparticle current j qp were considered in many papers starting from the Ginzburg paper [9] (for 
more references see [10]). These effects in a short bridge between two superconductors were studied recently in Ref. 
[10]. The authors established that due to a charge imbalance the thermoemf in the superconducting state may be 
comparable with that in the normal state. The thermoemf Vth induced in mesoscopic S-N structures with normal 
reservoirs kept at different temperatures was theoretically studied in Refs. [11-13,17]. Using a scattering matrix 
approach, Claughton and Lambert [11] studied the dependence of the thermoelectric voltage on the phase difference 
ip between the superconductors S taking into account the thermoelectric current jth = t]\7T. They showed that this 
dependence is periodic. 

Unlike Ref. [11] in Refs. [12,13] the thermoelectric current jth was completely neglected because the thermoelectric 
coefficient rj contains a small parameter: r) oc T/ep, where ep is the Fermi energy. It was shown that even in the absence 
of the thermoelectric current jth a temperature gradient in the N wire leads to the build-up of a thermoelectric voltage 
Vth if the N wire is in contact with two superconductors and the phase difference ip between the superconductors is not 
zero. The voltage V t h in the S-N structure may be of the order V t h ~ ci&Tje, where 5T = T r — Tj, T r j is the temperature 
of the right (left) normal reservoir (see Fig.l), c-i is a numerical factor (in the case under consideration ci ~ 0.1). 
Therefore the thermoemf Vth below T c is much larger than the ordinary thermoemf V or d (V or d ~ (T/ep)dT/e) above 
the critical temperature T c . This effect can be called a giant thermoelectric effect in S-N mesoscopic structures. It 
was shown in Refs. [12,13] that the voltage Vth oscillates with increasing the phase difference: Vth oc sin^, that is, the 
function Vth{tp) is shifted by 7r/2 with respect to the phase dependence of the resistance variation SRn of the normal 
wire: 5Rn(<p) oc cos p. The amplitude of these oscillations Vt/,(7r/2) is a non-monotonic function of temperature with 
a maximum at a temperature T m of the order of the Thoullcss energy e^. These results are obtained in Refs. [12,13] 
on the basis of microscopic equations for the quasiclassical Green's functions. It was assumed that the S-N interface 
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transmittance is low (due to a mismatch of the Fermi surfaces or to the presence of a potential barrier) and therefore 
the PE is weak. In this case the problem can be solved analytically. 

Another limit of perfect S-N interfaces was considered in Ref. [17] by using the same microscopic equations. Solutions 
for these (kinetic and Usadel) equations were found in Ref. [17] numerically and the obtained results are similar to 
those found in Rcfs. [12,13]. In the theoretical publications [12,13,17] the physics of this giant thermoelectric effect is 
explained in terms of the temperature-dependent Josephson current which, in the presence of a temperature gradient, 
has different values at different S-N interfaces. In order to compensate this difference, a quasiparticle current driven 
by the voltage V t h arises in the system. An approximate formula for the thermoemf was presented in Ref. [17], where 
in the main approximation the voltage Vth is expressed in terms of the Josephson critical current I C (T). A correction 
to this expression for the voltage V t h is small. Although in some limiting cases the representation of V t h through 
I C {T) is possible (the authors assumed that the order parameter A is much larger than the Thouless energy e L ), in 
a general case this can not be done. The point is that the Josephson current is a thermodynamical quantity (it can 
be presented as a derivative of a free energy with respect to the phase difference), whereas the voltage Vth is not. 
To make this point quite clear, one can consider a limiting case of temperatures T larger than cl (to be more exact, 
the ratio 2ttT/cl should be much larger than 1). In this case the Josephson current is exponentially small and the 
Josephson coupling is negligible. However, generally speaking, the voltage V t h is not exponentially small. Its value 
is determined by the ratio of other parameters: A and e^. The aim of this paper is to consider the case when the 
Josephson current is small, but the thermoelectric effect is not small in comparison with its maximal value. We will 
show that if the ratio 2irT/e L is large, the Josephson coupling between superconductors is almost negligible, but the 
thermoemf V t h is not small and can be measured. In this case one can say about a long-range thermoelectric effects. 



II. MODEL AND BASIC ASSUMPTIONS 



As in Refs. [12,13,17], we consider a system shown schematically in the inset of Fig. 3. A normal wire N, or a thin 
film with a width narrower than the coherence length £t = \J Dn/2ttT, connects two normal reservoirs N^ r . The left 

reservoir Nj is kept at a temperature T and the right reservoir N r has a temperature T + 5T. We assume that the S-N 
interface resistance R b is larger than the resistance Rl of the N wire so that the ratio i?i/i?f, is a small parameter. In 
this case the problem allows an analytic solution because there is a small parameter, the amplitude of the condensate 
function | F R ^ | , where F R ^ is the retarded (advanced) quasiclassical Green's function induced in the N wire due 
to the PE. Therefore the PE is assumed to be weak and the characteristics of the N wire deviate only slightly from 
those in the normal state. 

The thermoemf Vth arising in the system at ST ^ is an integral over energies e from a distribution function /-(e) 
[12,13,17]. This distribution function called a transverse part of the distribution function in Ref. [18] and denoted by 
fi in Ref. [19] is a difference between electron- and hole-like excitations: /_ = nj — px(scc, for example, [20] and [21]). 
This function determines the voltage and is related to a so-called charge-imbalance [22]. Another distribution function 
/+ = 1 — (raj +Pi) determines, for example, the supercurrent and the order parameter (in the superconductor). These 
distribution functions obey kinetic equations [18-21]. The kinetic equations were applied to the study of transport 
properties of mesoscopic S/N structures [23]. The function /_ obeys a kinetic equation, which being written in 
notations of Ref. [12,13] has the form 

M_ d x f-(x) +J s f+(x) - J an d x f+(x) = J (1) 

where all the coefficients are expressed in terms of the retarded (advanced) Green's functions: = G R ^a z + 

pR(A). M _ = (i_ g r G a - {F R F A ) 1 )/2; J an = {F R F A ) z /2, is an anomalous current, and the notation (F R F A ) Z 
means (F R F A ) Z = Tr{a z F R F A } /2. The function Js determines the Josephson current and can be expressed through 
the values of F R< - A ^ at the S/N interface 

J s = -(i/(4R b a))Tr{a z (F R F R - F A F A )} (2) 

The Green's functions Fg^ in the superconductors S are not disturbed by the PE due to a low S/N interface 
transmittance, and therefore in the right superconductor, they are equal to 

Fs (A) = cos(^/2) + ia x sinfa/2)] A/£f (A) (3) 
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where £ e = \/(e~± ij) 2 — A 2 ; a phenomenological parameter 7 describes a damping in the superconductors; the 
matrices in the left superconductor have the same form if one makes a replacement: ip —(p . The coefficient 

M_ is also expressed through the retarded (advanced) Green's functions F R ^ in the N wire, but the corrections 
due to the weak PE are small and approximately we have M_ « 1. The partial total current J is constant over each 
separate piece of the N wire. Since the "currents" Js and J an are proportional to small parameters | F R ( A ^ | 2 , the 
distribution function /+ (longitudinal in terms of Ref. [18]) may be taken in zero order approximation, that is, equal 
to its value in the absence of the PE. It has the form 

f + = [Sf eq /2}(l + x/L N ) (4) 

where Sf eq — [tanh(e/2(T + ST)) — tanh(e/2T)]. Eq.(l) should be solved with the boundary conditions at the points 
x = ±L N : f_(±L N ) = F_(V rt i), where F_(V) = [tanh((e + eV)/2T) - tanh((e - eV)/2T)}/2. In addition, one has 
to use the boundary conditions at the S/N interfaces (see [12,13] ). The voltages V± = (V r ± V{)/2 are found from 
the condition of the absence of the total current through the N reservoirs (the condition of a disconnected circuit). 
They are equal to [13,30] 

eV+ £* —5T(Li/L N ) [ de(e/3)g z+ (e, L x )fJ f de((v s + <?i+(e, L{))f' ); (5) 
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J de(e/3)g z _(e,L 1 )f' eq (6) 



where g z± = (1/4)[(F* T F A )(F R ± F*)] z , g 1+ = (l/4)[(i?« + F A )(F* + F A )] U f eq - cosh- 2 (e/2T), and v s = 
ReG^(e) is the density-of-states of a BCS superconductor, r s = Rs/Rb is a small parameter, R b is the S/N boundary 
resistance, Rs — L\/a is the resistance of a piece of the N wire in the normal state. The denominator in Eq.(5) is 
proportional to the conductance of the S/N interface [23-25]. The first term in the integrand of the denominator 
describes the conductance due to quasiparticles above the gap A, whereas the second term describes the subgap 
conductance (|e| < A). Therefore the integrand is not zero at all energies. 

Thus the voltages V± are expressed in terms of the Green's functions F R ^ which can be found from the linearized 
Usadel equation. Discuss now some general properties of the expressions (5) and (6). The condensate functions F R ^ 
induced in the N wire have the same matrix structure as the functions Fg^ (see Eq.(3)) 

F R = iby cos(cp/2)F R + ia x sm(cp/2)F R (7) 

Taking into account that F R = —(F A )*, the expressions for g z ± and for gi + can be represented in the form 

1 1 

9z+ = ^ Re(F y - F x ) Im F s sin <p; g z _ = - lm(F v - F x ) Re F s sin <p (8) 



5 i+ = - lm[(F y + F x ) + (F y ~ F x ) cos <p] Im F s (9) 



Here we dropped the index R: F R = F y . Note also that ImF s = -A/ | ^(e) | at | e |< A and Im F s = | at | e |> A, 
whereas ReFs = A/^(e) at | e |> A and ReFg = at | e |< A. Therefore at low temperatures (T << A) the 
voltage V- is much smaller than the voltage V+. It is useful to compare Eqs.(8) and the expression for Js that can 
be represented as 

J s = (r s /2Li) lm[(F v - F X )F S ] sin^ (10) 
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FIG. 1. Temperature dependence of the normalized thermoemf (the solid lines, the left scale) v+ = e(V r +Vi)/2ST and the nor- 
malized critical current (the dashed lines, the right scale) I C (T) / I c (0) for different values of the ratio A/es and rs = 0.5; where 
es = D/Li and A is the energy gap at zero temperature. Results are shown for L1/L2 = V%; Li/Lsn = Li/(Ljv — L\) = \/0.3; 
7 = O.les- 
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FIG. 2. The same dependencies as in Fig.l for 7 = 0.03es. 

The critical current is an integral over all energies from the function Js(e) which can be transformed into a sum 
over the Matsubara frequencies uj — nT(2n + 1) 

/c = S / deM{Fy ~ Fx)Fs] tanh ( £ / 2T ) = Et^ - F *) F s]<=i» (ii) 

We will see that at large ratio T/e^, the difference (F y —F x ) oc cxp(— 4y / 2a7/e7) is exponentially small and therefore 
the critical current I c is very small (here we take eh = D/L 2 with L = L\ = L2). On the other hand, Eq.(8) for g z+ , 
for example, can not be represented as a sum over the Matsubara frequencies. The point is that the critical current 
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I c can be written as an integral over energies from the product of the retarded (or advanced) Green's functions only: 
I c ~Im / de[(F y -F x )Fs] R tanh(e/2T). Therefore one can enclose the contour of integration in the upper half-plane of 
e, where the retarded functions are analytic functions, and calculate the sum over the poles (the Matsubara frequences) 
of the function tanh(e/2T). The functions g z ± contain the products of the type (F y — F x ) R Fg (an anomalous function 
in terminology of Ref. [26]), which are not analytic functions both in the upper and lower half-plane of e. Thus the 
integral can not be reduced to the sum over the Matsubara frequencies. The importance of the anomalous terms was 
emphasized by Gor'kov and Eliashberg [26] who showed that due to these terms the generalization of the Ginzburg- 
Landau equation to a nonstationary case, generally speaking, is not possible. In the case of mesoscopic S/N structures 
these terms lead to long-range effects [27]. 

Thus in a general case the statement about the smallness of V± is not valid. Only if the condition T, << A 
is fulfilled, one can regard F s as a constant (the integrand in Eq.(8) converges over energies of the order of T) and 
represent Eq.(6) as a sum over the Matsubara frequencies. In this case the voltages V± are also exponentially small. 
However if the Thouless energy is not small in comparison with A (for example, in [8] es ~ 2.8 K and A w 2.28 
K) , one can not neglect the dependence of Fs on the energy e and represent the integral in the form of the sum over 
the Matsubara frequencies. In this case the critical current I c may be exponentially small (J c oc exp(— y^rT/et)), 
whereas the voltages V± are not small. 

In order to calculate the voltages V± explicitly, we need to find the retarded (advanced) Green's function F R ( A \ 
In the considered limit of the weak PE, these functions are easily found from a solution for the linearized Usadel 
equation 

d 2 F R W/8x 2 - (^(^pRiA) = (12) 

where (n R ^) 2 — ^2ie/DM- We write out here the solutions for the functions F x y (we again dropped the indices 
R(A)) 

F x = {rsFs/0 y )[tanh9 s tanh9sN + tanh 0,, (tanh 5 + tanh0 SN )]/V x (13) 



F y = (r s Fs/0y)[tanh6sN + tanh0„(l + tanh0 s tanhflsjv)]/^ (14) 

where V x = tanh0 a tanh0s tanhOsN + tanh0s + tanhOsN', 

V y = l + tanh0sAr(tanh0 y + tanh0s); r s = Ri/Rb, Ri = Li/a, y = nL 2 , 9 S = kLi, 9 S n = kL S n, L S n = L N — Li. 

The exponential smallness of the critical current I c at high temperatures T > D/L 2 can be easily verified if one 
considers a simple case of equal distances L\ = L2 = Lsn = L, that is, 9s = y = 9sn = 9. In this case we get for 
the difference (F y - F x ) 

F.-F.- (rsFsm^ i2 + t ^l^ e) (15) 

It is seen from Eq.(15) that at >> 1 we have (F y - F x ) oc (1 - tanh 2 0) 2 cx (1/0) exp(-4Ly / 27rT/D), that is, this 
difference is small if 2co w 2irT >> D/(2L) 2 . Therefore the critical Josephson current also is exponentially small (see 
Eq.(ll)). On the other hand, neither V+ nor V- can be represented as a sum over the Matsubara frequencies. Thus, 
one can not claim that, for example, V + should be exponentially small if the Josephson coupling is very weak. 



III. RESULTS AND DISCUSSION 



We can rewrite the expressions for V± in the form 

eV+/ST = -(LJLn) < {e(3)F'_Fg > e sin^/ < 2v s + [F" + F" cos V }F s ' > £ ; (16) 

2eV-/ST = r s < F'iF' s > e sin<^. (17) 

where rs = Li/R^cr, F' ± — Ke(F y ± F x ), F' s ' — ImFs etc. The angle brackets mean the averaging over energies: 
< ... >= / de{...)f , and the functions F x ^ y are given by Eqs.(13,14). Using these formulae, one can numerically 
calculate the thermoemf as a function of temperature for different values of parameters. 
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In Fig.l and Fig. 2 we plot the dependence of the normalized voltage v + = eV+/8T and the normalized critical 
current i c = Z c (T)// c (0) as a function of the normalized temperature T/es- These curves are presented for (p = n/2, 
where v+ reaches a maximum, and different values of the ratio A/es. We choose the value of rs equal to: rg = 0.5 
(note that the curves for smaller rs, rs < 0.3, are similar to those plotted in Fig.l and 2). We set 7 = O.les for 
the data in Fig.l and 7 = 0.03es for the data in Fig. 2. Note that the damping parameter 7 may essentially exceed 
the corresponding value in bulk superconductors because the PE leads to an additional damping of the order of 
arses, where a is a parameter which depends on geometry of the S and N films. We see that in both figures the 
critical current decays to zero smoothly and fast (faster than exponentially) when the temperature T increases, but 
the magnitude of v + decreases not so fast. Close to the critical temperature T c the critical current is very small in 
comparison with J c (0), whereas the normalized thermoemf is not so small. 




FIG. 3. Temperature dependence of the normalized thermoelectric voltage between the normal reservoirs v- = e(V r — Vi)/25T 
for different values of the ratio A/es and rs = 0.5 (the values of other parameters are the same as in Fig. 2). Inset: the 
structure under consideration; the normal wire connects two normal reservoirs kept at temperatures T (left reservoir) and 
T + 5T (right reservoir) and two superconductors with the phases ip/2 and —<fi/2. The distance between the normal reservoirs 
is: Ljv — Li + Lsn- The electric potential of the superconductors is set to zero. 



Although the value of v_ = eV-/5T is much smaller than v + , it depends on temperature quite differently being 
very small at low T and having a maximum at a temperature near T c . We represent the temperature dependence 
of v_ in Fig. 3. As we mentioned before, the voltage eV- is very small at low temperatures T because the function 
g z -(c, Li) is zero at |e| < A and grows for T > A. At T = T c this voltage as well as V+ turns to zero. Therefore, for 
some parameters the voltages V r> i — (V+ ± V—) may have two maxima (or extrema): one at low T, where V+ has a 
maximum, and another one at a temperature close to T c , where V— has a maximum. It is worth mentioning that a 
possible assymetry of the system (for example, different distances between the crossing point and the left and right 
N reservoirs) may increase the voltage V- [17]. 

In Fig. 4 and 5 we plot the temperature dependence v r = v + +u_ for the damping parameters 7 = 0.1 and 7 = 0.03. 
One can see that this dependence has two maxima. The resistance of a S/N system has a similar temperature 
dependence with two maxima [28]. Moreover, for A/es 1-5 the thermoemf drops to zero at a certain temperature, 
remains very small at larger T and increases again at temperatures close to T c . At T = T c the thermoemf drops to 
zero. A similar dependence of the thermoemf V r as a function of the magnitude of the dc heater current was observed 
in experiment [8]. 

Note that the voltage Vi = V+ — V- of the left cool N reservoir (see Fig. 6) has the second maximum near T c the 
sign of which is opposite to the sign of the first maximum . Such a behaviour was observed in experiment [29] . 

Obviously the average temperature in the N wire (or film) is proportional to the magnitude of the heater current. 
Thus one can say about a qualitative agreement between theory and experiment. It is difficult to carry out more 
precise comparison because the PE in the experiment seems to be strong. Perhaps, the strong PE is the reason for 
a high value of the second maximum of V r near T c (in the experiment, it is comparable with the value of the first 
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maximum at low T) because it is proportional to the parameter rs- 




0.14 r- 
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FIG. 5. The same dependencies as in Fig.4 for 7 = 0.03es. 

We assumed that the parameter rs is small; however in the case of a strong PE this parameter should be taken of 
the order of 1 because the condensate functions F R ^ in the N film saturate with increasing rs- Finally we give some 
estimations for parameters close to the experimental ones [8]: for L\ — and D — 150cm 2 /s, one has eg w 2.&K. 
This value of es is a little larger than A « 2.2K (A « 0.865). The estimations of the thermoemf give the values 
eVth/ST w 6 — 10/jV/K. These values are comparable with the voltage V t h observed in Ref. [5], but much larger than 
Vth measured in Refs. [6,8] (note that in Refs. [5,6] the voltage V-, but not V+ was measured; in the case rs ~ 1 
the magnitudes of both voltages are comparable). Perhaps the reason for this discrepancy is the energy relaxation 
processes in the N wire which were neglected in our calculations. Indeed, in order to neglect the inelastic scattering 
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the condition r^ 1 << rs(D/L 2 ) should be fulfilled, where t^ 1 is the inelastic scattering rate. Otherwise in the system 
arises a strong depairing which suppresses the PE and destroys the phase coherent effects discussed above. The more 
exact comparision with experimental data is not possible at this stage; for example, a thermoemf symmetric with 
respect to the phase difference ip was observed in experiments. The origin of this part of thermoemf stll is unclear. 
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FIG. 6. Temperature dependence of the thermoelectric voltage vi — (v+ — V-)= eV r /ST for different ratio A/es , for rs = 0.5 
(the values of other parameters are the same as in Fig. 2). 



IV. CONCLUSION 



In conclusion, we show that the thermoemf arising in the four-terminal S/N structure (see inset in Fig. 2) is measur- 
able even in the case of a very weak Josephson coupling. One can say about a long-range thermoemf. The thermoemf 
Vi tr (or thermopower Vi ir /ST) depends on T in a non-monotonic way and may have two extrema: one at low tem- 
peratures and another at temperatures close to T c . At some values of parameters and intermediate temperatures the 
thermoelectric voltage is negligible. This behavior qualitatively agrees with the experimental observations [8] . 
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